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Abstract 

Gamma ray earthbound and satellite experiments have discovered, over the 
last years, many galactic and extra-galactic gamma ray sources. The de- 
tection of astrophysical neutrinos emitted by the same sources would imply 
that these astrophysical objects are charged cosmic ray accelerators and help 
to resolve the enigma of the origin of cosmic rays. A very large volume 
neutrino telescope might be able to detect these potential neutrino emitters. 
The apriori known direction of the neutrino source can be used to effectively 
suppress the 40 K optical background and increase significantly the tracking 
efficiency through causality filters. We report on advancing filtering and pre- 
fit techniques using the known neutrino source direction and first results are 
presented. 
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1. Introduction 

The future Mediterranean very large volume neutrino telescope, KM3NeT, 
will be one of the world's largest particle detectors and will provide a research 
infrastructure for a rich and diverse deep-sea scientific program [H,0]. Various 
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astrophysical sources are expected to produce high-energy neutrinos that can 
be detected with KM3NeT. The most promising candidate neutrino sources 
are Galactic point-like (or with a small angular size) sources. The obser- 
vation of neutrinos emanating from such sources would bring unique new 
insights on the nature of cosmic accelerators and might help to resolve the 
enigma of the origin of cosmic rays. Observations by gamma ray earthbound 
and satellite experiments have revealed many such astrophysical objects, in 
which high-energy processes at and beyond the TeV energy level take place. 
However, measurements with gamma rays alone cannot clearly distinguish 
whether the accelerated particles are leptons or hadrons. Only the obser- 
vation of neutrinos from a source can unambiguously establish the hadronic 
nature of that source. 

However, the Galactic sources are generally expected to have a cut-off in 
their energy spectra in the range 1 — 10 TeV. The signal produced by neutrinos 
in this energy range could be lost among the 40 K optical background, result- 
ing in low track reconstruction efficiency In this work we report on advancing 
triggering, filtering and prefit techniques that use the apriori known direction 
of the neutrino source to effectively suppress the 40 K optical background and 
enhance the detection efficiency for neutrinos from Galactic sources for which 
the direction is known from their gamma ray emission. 

2. Detector description and simulation framework 

In the present study the telescope layout considered is the one optimized 
during the KM3NeT Design Study jij, exhibiting optimal sensitivity in dis- 
covering astrophysical point sources emitting neutrinos with an energy spec- 
trum of E~ 2 and either a high energy cut-off or no cut-off at all. According 
to this layout, the KM3NeT detector will consist of 12320 photo-sensors dis- 
tributed over 308 Detection Units (DUs). The photo-sensor unit is a digital 
optical module (DOM) consisting of a 17-inch diameter pressure resistant 
glass sphere housing 31 3-inch photomultiplier (PMT) tubes jijj]. A DU is 
a vertical structure which carries photo-sensors and devices for calibration 
and environmental measurements, arranged vertically on 20 Storeys. Each 
Storey consists of a bar with one DOM at either end. The horizontal dis- 
tance between the centers of the DOMs is 6 m. The vertical distance between 
Storeys is 40 m, while the position of the lowest Storey is 100 m above the 
seabed. The bars have an orientation orthogonal to their neighbors. The 
distribution of the positions of the DUs on the seabed (the so-called foot- 
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print) is homogeneous. The footprint forms a roughly circular shape and has 
a typical DU density corresponding to an average distance between neigh- 
boring detection units of about 180 m. The total instrumented volume of 
the detector is 5.8 km 3 . 

For the present study we have simulated the detector signal produced by 
neutrinos using the HOURS (Hellenic Open University Reconstruction and 
Simulation) physics analysis package Q. We have used a generic neutrino 
flux distributed isotropically (in a 4-7T solid angle) on the Earth's atmosphere, 
with an energy distribution following a power law spectrum, with a spectral 
index of —2.0, in the range of 15 GeV — 100 PeV. 

In the following, the signal produced by the secondaries of a neutrino 
interaction event in the vicinity of the neutrino telescope is described by a 

— * 

collection of active DOMs (hits), hi(ti, H^mi) (i = l,..,Nhu)- Each hit is 
described by the first photon arrival time, the active DOM's position, H^, 
and the multiplicity, m,. The hit multiplicity is the number of active 3-inch 
PMTs of a DOM within a time window of 10 ns. A level one (LI) coincidence 
trigger is defined to fire/be activated for rrii > 1. 

A random background rate of 5 kHz is assumed for each 3-inch PMT, 
including dark current, 40 K decays, and bioluminescence. In addition to 
random coincidences, an LI hit (with multiplicity m = 2) rate of 500 Hz is 
assumed on each DOM, due to genuine coincidences from K decays. 

For each signal event, produced by a neutrino emanating from a known 
astrophysical point source, we assume that the direction vector, d, of the 
muon track coincides with the source direction. In order to estimate the fake 
signal produced by implying in the analysis a certain direction, we assume for 
background events (atmospheric neutrinos and atmospheric muon bundles) 
a random (fake) candidate muon track direction distributed isotropically in 
a 4-7T solid angle. 



3. Background filtering 

The filtering technique of the m K optical background is based on a causal- 
ity criterion that uses the assumed direction vector, d, of the muon track 
generated by the incident neutrino. The arrival time, ti, of a photon emitted 
by the muon with the Cherenkov angle, 9 C , to a DOM located at position Hi 
satisfies the relation: 

ctj = Qj + bitan9 c , (1) 
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where dj = d- (Hi — V), and bi = \Hi — V — did\ is the vertical distance of the 
DOM to the muon track. V is the position vector of the pseudo-vertex, cor- 
responding to the position of the muon when the time measurement started. 

Two hits hi(ti, Hi, mi) and hj(tj, Hj,mj) produced by direct photons on 
two different DOMs according to Eq. [TJsatisfy the following relation, in which 
V is eliminated: 

cAt -d-AH 

n = Ab , 2 

tan9 c K ' 

where At = U — tj, AH = Hi — Hj and Ab = hi — bj. If we project the 
DOM's positions on a plane perpendicular to the assumed direction d, then 
from simple geometry and using Eq. [2] we have the relation: 



IA6I 



cAt-d- AH 



tanOr 



< \AH - {d ■ AH)d\. (3) 



According to Eq. [3] we form the causality criterion to be satisfied by the two 
hits: 

\cAt - d ■ AH\ < tan9 c \AH - (d ■ AH)d\ + ct s , (4) 

where t s = 10 ns is inserted in order to relax the causality requirement, due 
to light dispersion and time jitter effects. In order to limit the number of 
noise hits satisfying this criterion, we also require the longitudinal distance 
between the two DOMs along the direction of the muon track to be: 

\d- AH\ < 800 m, (5) 

while the corresponding lateral distance should be less than one absorption 
length^: 

\AH - (d ■ AH)d\ < 67.5 m. (6) 

Equations IU [5] and |6] define the causality condition that has to be satisfied 
by any two hits. 

In order to use this condition as criterion for background filtering, for 
every hit, hi, we calculate the following sum: 

N hii 

Mi = J2 m 3 C ( h h h j) 



This is the maximum absorption length for wavelength 440 nm. 
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where rrij is the multiplicity of the hit hj defined in Section HI while C(hi,hj) = 
1 if the causality condition is satisfied by the two hits, or otherwise. If 
Mi < 5 the hit hi is rejected. When we apply this filtering technique we 
achieve a 99.7% rejection of the noise hits, while more than 90% of the signal 
hits survive. 



4. Triggering and prefit 



For every three hits on different DOMs, that when combined in pairs 
satisfy the causality condition, a pseudo- vertex, V, can be found analytically, 
as described in the following. Using Eq. [TJwe calculate that: 



cti = d- (Hi - V) + bitan6 c 

d + z ' = bi (i = 1..3), 



where 



and Cj 



d-V 

tanO r 



(7) 



(8) 



cti — d ■ Hi 
tanO c 



Then, we define a new coordinate system with the z-axis parallel to the 
assumed direction d, the origin at the point defined by Hi, and the xz-p\&ne 
defined by the z-axis and the point defined by H 2 . The base vectors of this 
coordinate system are: 



x 



-(H 2 -Hi) 
V 

-d x (H 2 - 
P 

d 



1 

p L 



d-(H 2 -Hi) 



d 



Hi) 



In this frame, the projections of the three hits that we consider on the xy- 
plane are: Hi(0, 0), H 2 (p, 0), H 3 (q, r), where the values of p, q and r can be 
estimated easily from the positions of the hits. The unknown pseudo-vertex 
V(x, y, z) is expressed as follows: 

V = Hi + xx + yy + zz (9) 
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where from Eq. [HJ 

z — d-(V— Hi) = z'tan6 c -d-H x (10) 

Taking into account that hi in Eq. [7] is the vertical distance of the hit position 
from the track, then the x, y, z' satisfy the following set of equations: 

b l = {ci + z 'f = x 2 + y i 

bl = {c 2 + z'f = {x-pf + y 2 

b 2 3 = (c 3 + z') 2 = { X - q f + {y-rf 

with p, r ^ 0. The above set of equations and Eq. [10] can be used to calculate 
the explicit values of x, y, z, and then from Eq. [9] the pseudo- vertex. 

The probability to find three LI noise hits^l producing a pseudo- vertex has 
been found in our simulation studies to be of the order of 10~ 5 . Consequently, 
the requirement to find at least one pseudo-vertex using only the LI hits of 
an event, offers a good directional selection criterion (trigger) when searching 
for muons arriving from a certain direction. 

Obviously, in a signal event, the above procedure finds many candidate 
pseudo- vertexes using different hit triplets. Then we examine the distribution 
of these pseudo- vertexes in space searching for clusters. The cluster with the 
most participants, within a radius of 5 m, is chosen to estimate the position 
of the track pseudo- vertex. This track pseudo- vertex is defined as the average 
position of the participants in the cluster. We have found, using the simulated 
neutrino events described in Section [2J that this prefit technique estimates 
the true track pseudo-vertex with an accuracy of 2 m, as it is presented in 

Fig.m 

The estimated track pseudo-vertex and the assumed direction, d, can be 
used in order to reduce further the contribution of noise hits. For this we 
require the time residual of the hits to be less than 150 ns tf p \ < 150 ns), 
where the expected time can be derived according to Eq. [Uwhen using for V 
the position of the estimated track pseudo-vertex. Furthermore, we require 
that the distance of the selected hits from the estimated track (expressed by 
the vector d and the track pseudo-vertex) to be less than 4 light absorption 
lengths in the sea water (~ 270 m). After this track prefit stage the 99.97% 
of the noise hits have been rejected, while 90% of the signal hits remain. 



2 Produced by PMT dark current, 40 K decays and bioluminescence. In Section [6] noise 
from atmospheric neutrinos and muons is examined. 
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Figure 1: The cumulative distribution of the distance between the estimated pseudo- vertex 
and the true track pseudo- vertex in meters (solid line) . The median of the distribution is 
2 m. The corresponding distribution for the alternative prefit technique used in HOURS, 
where an apriori known track direction is not used, is also shown (dashed line). 

5. Muon track reconstruction and results 

We use a combination of a \ 2 fit and a Kalman Filter selection algorithm 
in order to reconstruct muon tracks as described in j^j. The hits in the track 
reconstruction are those surviving the prefit stages described in Sections [3] 
and HI At this tracking stage the assumed direction d and the estimated 
track pseudo- vertex are used as initial values for the fit, however the track 
parameters are not constrained to the assumed direction. 

This tracking algorithm produces many candidate tracks. We select the 
best candidate track by employing information offered by the hit arrival times 
and Multi-PMT direction likelihood as described in (^j. Specifically, for each 
candidate track we evaluate the likelihood value, using the hits that have 
survived the background filtering stages. The best candidate track, with the 
largest likelihood value, is further improved by performing a final likelihood 
fit. However, the estimated track is finally accepted when it is misaligned 
with respect to the assumed direction vector, d, by less than 1^- 



3 For point source searches this angular difference has to be further optimized. 
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Figure 2: Reconstruction efficiency versus neutrino energy for events with at least three 
signal LI hits (solid line). The corresponding efficiency for the alternative fitting technique 
used in HOURS is also shown (dashed line). The lower values of the solid line than those 
of the dashed line at 7.5 are a statistics effect. 



for simulated events, containing at least three signal LI hits, produced as 
described in Section [2j The corresponding efficiency of the alternative fitting 
procedure used in HOURS (when the apriori known track direction is not 
used at all) is also shown. The reconstruction efficiency is improved by a 
factor of 2.5 for 1 TeV neutrino induced events, and a factor of 1.5 for 10 TeV 
neutrino induced events. 

Fig. [3] presents the point spread function for reconstructed events, con- 
taining at least three signal LI hits in comparison with the alternative fitting 
technique used in HOURS. The point spread function is the cumulative dis- 
tribution function of the angular difference between the incident neutrino 
direction and the reconstructed muon direction. The two techniques, for 
E~ 2 neutrino spectrum, exhibit the same resolution as is presented in Fig. [3] 
(the point spread functions have the same median of about 0.2°). 



The tracking quality criteria applied are minimal. We only require more than seven 
hits to participate in the best candidate solution found by the tracking algorithm. 



In Fig. |2]the reconstructio: 




versus the neutrino energy is shown 
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Figure 3: Point spread function for the reconstructed events with at least three signal 
LI hits (solid line). The corresponding point spread function for the alternative fitting 
technique used in HOURS is also shown (dashed line). The y-axis shows the fraction of 
events with an angular difference less than the given value. 

6. Estimation of fake signal 

In principle, taking into account the known direction of the source could 
produce fake signals. In order to estimate the probability of such a fake signal 
production we have simulated background events (atmospheric neutrinos and 
atmospheric muon bundles) assuming to be produced by a hypothetical point 
neutrino source with a random position isotropically distributed in a 4n solid 
angle. We applied on the signal produced by these background tracks the 
filtering, prefit and track reconstruction procedures described in Sections EJ 

Hand ED 

Fig. H] presents the probability of an atmospheric neutrino induced event 
to be reconstructed so that the angle between the reconstructed track and 
the assumed direction of the source is less than 1°. This probability is shown 
as a function of the cosine of the angular difference, 8, between the true at- 
mospheric neutrino direction and the assumed direction of the source. The 
probability to produce a fake signal is negligible for cos(8) < 0.9, i.e. atmo- 
spheric neutrino background has less than 10~ 3 probability to produce fake 
signal when the incident atmospheric neutrino comes from a direction which 
differs more than 25° from the point source. 
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Figure 4: Probability of an atmospheric neutrino induced event to produce fake signal 
versus the cosine of the angular difference, 9, between the true neutrino direction and the 
assumed direction of the source. 

However, this fake signal can be further reduced by applying adequate 
tracking quality criteria. Fake tracks have a low tracking quality which is 
reflected on the estimated tracking error, i.e. fake tracks carry a very small 
weight in the fit technique pj. 

7. Conclusions 

We developed a triggering, filtering and prefit algorithm taking into ac- 
count the apriori known direction of the neutrino source. The performance 
of the KM 3NeT detector in detecting/ discovering point-like neutrino sources 
using this technique is under further evaluation. The performance of the 
combination of the known source direction technique with the fit technique 
enhances the discovery potential of the detector (a publication is under prepa- 
ration). 
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